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AESTRACT 



The Generalized Rotor Performance (GRP) program is a 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used tc evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. The program was originally 
designed by the Sikorsky Aircraft Company and purchased by 
the United States Navy. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to add selected 
desirable features to the program, (3) to produce a much 
needed Users' Manual, and (4) to run an analysis comparing 
the program's calculated results against manufacturer's 
data. These goals were accomplished and the results of the 
analysis indicated that the program produces highly accurate 
results within the normal cruise range of a modern 
he licopter . 
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INTRODUCTION 



The Generalized Rotor Performance (GRP) program is a 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used to evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. The GRP series of programs 
were originally designed by the Sikorsky Aircraft 
Corporation with the cooperation of the United Aircraft 
Research Laboratories. Work began on this series of 
programs in 1955 and has continued to date. The United 
States Navy purchased this program in 1964 and has used it 
in the Naval Air Systems Command as a helicopter performance 
computer program. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to ensure that the 
Navy's version of this program performed the calculations 
according to the design theory, (3) to produce a much needed 
Users' Manual, (4) to add certain desirable features to and 
correct errors found in the program, (5) to run an analysis 
comparing the program's calculated results against 
manufacturer's data and (6) to document areas that may need 
further attention. The work was done with the help and 
cooperation of both faculty at the Naval Postgraduate School 
and personnel in the Naval Air System Command. 

The output available from the program includes rctor 
shaft horsepower, rotor profile horsepower, main rotor 
torque, lift, drag, thrust, H force, rolling and pitching 
moments and other forces and moments calculated in the 
shaft, control and relative wind reference axis systems. 
Also available are the azimuthal histories of the blade's 
flapping angles, rates and accelerations, plus azimuthal and 
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radial histories of angles of attack, CL, CD, Mach number, 
sweep angles and air load distribution. In addition, the 
program will calculate a Fourier coefficient series for 
flapping angle, radial station air leads and azimuthal Z 
force distribution. 

The program can be divided figuratively into three main 
routines. The first routine determines a steady state 
flapping solution. The second routine determines the forces 
and moments generated by this flapping solution. The third 
routine compares the calculated results of lift, drag and 
other options to those desired by the user. If the 
calculated values are not within a predetermined tolerance, 
the program will, by use of a first order Taylor series, 
generate new values to reenter into the blade flapping 
routine. 

The program uses a combined blade element/momentum strip 
theory. The blade's radius is divided into a maximum of 15 
segments. The blade is advanced around the azimuth in a 
prescribed number of degrees. At each particular azimuthal 
position, velocities are computed in the perpendicular (UP) , 
tangential (UT) and radial (UR) directions. From these 
velocities and a known pitch angle distribution, the local 
angles of attack are calculated. From these angles of 
attack, the forces at each radial blade segment are 
determined. Once the forces are known, the local moments 
can be found. From the summaxion of these moments at a 
particular azimuthal position, the blade flapping 
acceleration is calculated. This acceleration, combined 
with the calculated flapping angle and rate, is used to 
advance the blade to the next azimuthal position. A steady 
state flapping condition is assumed tc exist when the blade 
flapping angle and rare at the PSI egual zero and 360 degree 
azimuthal position are within a prescribed tolerance of each 
other. 
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This method of calculation accounts for retreading blade 
stall, the reverse flow region and compressibi«iit v effects. 
The program can accommodate a full range of geometric and 
design variables including flapping hinge offset, elastic 
flapping hinge restraint, first and second harmonic cyclic 
inputs and spanwise variations in blade twist, local mass 
densities, chord and tip sweep. The program uses no small 
angle assumptions and has no restrictions on tip speed, 
forward velocity or advance ratio. The rotor system can be 
oriented in any direction in space and can be given any 
rotor shaft or aircraft roll, pitch or yaw angular 

velocities. The program will perform calculations in 

straight and level flight or a uniform induced velocity may 
be added to simulate climbs and descents. The GRP will 

accept up to five different airfoil data decks plus one 

blade spar characteristics data deck. Lift and drag 

information is entered into the program in the form of CL 
and CD versus Angle of Attack Tables for up to 15 different 

Mach numbers. The user has a choice of six methods of 

solution depending upon the desired restrictions. The user 
has available nine printout options, two of which are 
program debugging options. There are also error printout 
messages that will assist the user having difficulty with 
the program. 

In the method described above, the flapping motion 
determined is about a flapping hinge offset, and with the 
elimination of all assumptions that the flapping angles are 
small, the rotor control axis (axis of no feathering) and 
the tip path plane axis (axis of no first harmonic flapping) 
are no longer considered convenient for reference in the 
analysis of rotor blade motion. Therefore, the axis 

selected for use in the analysis is the rotor shaft axis 
system. All forces, moments and angles are referred to the 
shaft axis system with the exception of some of the final 
output forces which are referred to the relative wind and 
control axis systems. 
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With the use of the computer, the GHP program eliminated 
many of the simpliflying assumptions of the earlier 
classical theory originated by Wheatley and Bailey in Ref. 1 
and 2. The assumptions of the classical theory did not 
impose serious limitations in low speed flight, but as 
helicopter speeds increase, the inaccuracies inherent in the 
theory seriously limit the usefulness of this method. The 
assumptions of the Wheatley-Bailey theory which are not 
present in the GRP program include the following. 

1. The flapping and inflow angles are assumed 
small . 

2. The lift and drag coefficients are approximated 
by a linear and a quadratic variation, 
respectively, with angle of attack. 

3. The effects of Mach number on CL and CD are not 
considered . 

4. Sectional characteristics in the reverse flow 
region are the same as those in the 
conventional flow. 

5. The sectional characteristics of blade twist, 
tip sweep, flapping hinge offset and root cut 
out are ignored. 

While the GRP program represents a refined approach to 
the classical theory, there are still some basic assumptions 
which make the GRP subject to error. These are: 

1. Steady state two-dimensional airfoil data are 
used . 

2. Quasi-static blade analysis is used. 

3. The rotational speed about the shaft axis is 
constant. There is no lead cr lag motion in 
the program. 

4. The rotor blade is assumed rigid in bending and 
torsion . 



12 



“JUL 

SAV 



5. 



The rotor inflow is assumed uniform unless 
varied by the user. 

6. Spanwise flow is incorporated into the 
calculation of blade angle of attack; however, 
it is assumed that the flow at one segment does 
not affect the flow at any other segment. 

It is felt that the errors induced by the above 
assumptions are relatively small in the normal cruise speed 
region of a modern helicopter. This is the region from just 
below the minimum power airspeed to the maximum allowable 
cruise speed. The program can- not calculate rotor hover 

power. This is because one of the factors in the 
denominator in the main routine which estimates new flapping 
routine reentry parameters is the advance ratio. Since the 
advance ratio is zero in a hover, the computer would attempt 
to divide here by the number zero. Highly accurate results 
are not available in the airspeed region from hover to just 
below the minimum power airspeed while using the normal 
uniform inflow assumption. This is a region of highly 
non-uniform flew. Variable inflow may be inputea by the 
user. However, at this time, no information is available on 
how successful this technique is in accurately estimating 
the required rotor system horsepower. While the program's 
required technique for inducing harmonic variable flow is 
cumbersome, Bramwell in Chapter Four of Ref. 6 describes 
several methods which could be incorporated into the 
program. 

The quasi-static blade analysis assumption is valid 
except in the very high speed region where a large 
percentage of the retreating blade is in a stalled 
condition. This is an area of aerodynamic hysteresis that 
is influenced by unsteady aerodynamics. While encountering 
lift hysteresis, the amount of lift change above the steady 
state CL is about the same as below. However, the lift 
distribution on the rotor disc will change. Stall is 
delayed and occurs at a slightly later azimuthal station 
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then predicted by quasi-static analysis. While difficult, 
complicated and computer-time-consuming procedures can be 
taken to reduce these assumptions there is no guarantee that 
the accuracy of the solution will improve. It is felt that 
the present program produces highly accurate results in the 
flight range of interest in a modern helicopter. 
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METHOD OF AN ALY SIS 



A. GENERAL COMMENTS 

The GRP uses a combined blade element/momen turn strip 
theory in its calcultions. The blade radius is divided into 
a maximum of 15 segments or strips. Figure 1 illustrates a 
typical blade element. In order to obtain a steady state 
flapping solution the program must calculate the time 
history of the flapping motion. This necessitates a 
complete knowledge of the flapping angles, rates and 
accelerations. The time history cf this motion must be 
found by solving the highly non-linear equation for the 
summation of moments about the flapping hinge. Figure 2 
illustrates the forces which produce moments on the rotor 
system. The blade has aerodynamic forces acting upwards, 
blade weight acting downwards, centrifugal force acting 
outwards, an optional elastic hinge restraint acting to 
return the blade to a zero flapping angle and an inertia 
force resisting any change in blade flapping. The moment 
summation equation about the flapping hinge may be written 
as 

(') * "» "s * « tf * 

• • 

The moment due to inertia force can be expressed as 

where I t is the blade moment of inertia about the flapping 

• • 

hinge. Substituting this into equation 1 and solving for 
the following governing equation is obtained. 

(2) ^ 

If the angle of attack and local Mach number at each 
blade segment were known, the CL and CD could be obtained 
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Figure 1 Blade Element Diagram 
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from the inputed airfoil information. The local aerodynamic 
forces and moments can be calculated and the flapping 
acceleration determined. Equation 2 can be rewritten as 



(3) $ 



n 



I ( 

i-\ 



+ dM^ + ) 



The flapping angle and rate at the (II + 1) azimuthal 

position could be obtained from the following expressions. 



( 4 ) ?rv^v - ^ ^ 

X 

(5) 4 . '<j A 

However, since the flapping is periodic in nature and 
has a direct relationship to the azimuthal angle, PSI, the 
values for flapping are solved with respects to PSI, vice 
time. The values of SI, y and time are related by the 
equation =-A&t. Therefore 



(6) 

(7) 



• • 



=. SlA& tr. -H-& 

at. dv 
£± = 

dt? - 6 ^ 



The governing equation for the flapping motion now 
becomes 



1 %. 

( 8 ) § = 



The flapping solution is based on the assumption that 
the angle of attack is known. However, it is not and the 
program must proceed through an iterative process in order 
to determine the inflow ratio, collective pitch and cyclic 
input angles required to generate the desired forces. 
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ANGLE OF ATTACK CALCULATIONS 



A very important and basic part o 
procedure by which the local an 
calculated. While the pregram will 
attack with any angular velocity appli 
system or the helicopter, the developm 
level flight only. The classical ap 
flow, UR, and the angle of attack woul 
in Figure 1. However, as the blade ro 
it will encounter a large variation 
program attempts to compensate for th 
following manner. Instead of the infl 
the arctangent of UP/UT, it is set e 
UTUR is the resultant velocity in the 
direction. This is illustrated in 
angle is also reduced by the cosine cf 
angle of attack is now calculated in 
three-dimensional angle of attack 
classical two-dimensional angle. 

The program enters the CL, CD t 
plane angle of attack and the sweep 
number. The program computes the fore 
in the sweep plane, UP and UTUR, 
geometry in the normal plane. Once t 
in the sweep plane they are resolved 
directional forces. 

This three-dimensional angle of 
will delay stall on the rctor by r 
attack. This describes what actually 
However, it is felt by previous perso 
that there is a point where the sweep 
it tends to wash out the lift being 
has been modified to reduce CL b.y one 
between 60 to 72.5 degrees and reduce 
angles between 72.5 and 90 degrees, 
occur normally only on the inboard 



lft. 



f this program is the 
gles of attack are 
calculate angle of 
ed either to the rotor 
ent here will describe 
proach ignores radial 
d be obtained as shewn 
tates about the shaft, 
in radial flow. The 
is radial flow in the 
cw angle PHI equalling 
qual to UP/UTUR where 
tangential and radial 
Figure 3. The pitch 
the sweep angle. The 
the sweep plane. This 
is lower than the 

ables with this sweep 
plane resultant Mach 
es using the velocites 
and the blade chord 
he forces are computed 
into their respective 

attack, due to sweep, 
educing the angle of 
occurs on the blade. 
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Spanwise Flew Diagram 
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retreating blade in and near the reverse flow region. 

DP, UT, UR and pitch angle for level flight are shown in 
Figures 4 and 5. They are calculated in the following 
manner. The radial velocity, UR, is calculated as 



(9) UR = ( Vcos^ ) cos H' cos ^ 



The tangential velocity, 
first is the local rotational 
sinusodial component of the 
general expression for UT is 



UT, 


has 


two 


components. 


The 


velocity 




the second 


is a 


forward 


flight 


velocity . 


The 



(10) UT = Jt r + (Vcoso^j) sin^ 

This expression contains a small angle assumption in the 
term jir for blade flapping angle. The program accounts for 
the fact that the flapping angle does reduce the true radius 
slightly by using the following formula 

(11) UT = (E/R + (r - E/R) cos ^ )/L + (Vcos^) sin^ 

The perpendicular velocity, UP, consists of three terms, 
the inflow ratio, a flapping velocity and a small component 
cf forward velocity. The inflow ratio is defined as 



(12) 


= 


VsinVa 


- Or 
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The 


second 


component is a v 


eri 


cal flapping 


velocity 


which is 


a function of 


flapping r 


ate 


and 


radius. 


This 


is 


computed 


as 
















(13) 
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= (r - 
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UT = (E/R +• (r - E/R)cos$Ux + (Vccs 0:5)3 in 4> 
Figure 4 LfT Diagram 



21. 



CT 








UP - A/lRcos$ + (r - E/R)(3 + ( Vcosdi )cosYsin(^ 



UR - ( Vcosd* ) cos'f cos@ 






(b) 




(.Vco-s^ s'i co^y 




Figure 5 UP and UR Diagram 
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(14) FJP (3) = (7cos«( s )cos1' sin?) 



The total formula for OP is 



( 15 ) 



UP = )i4i.Rcos^ + (r - E/H)(J> + (Vcos^s) cosfsin^ 



The pitch angle (0) is expressed in equation 16. 0.75 
is the pitch angle at the 75 percent radius position, 0 / is 
the twist depending on the relationship of the location of 
the blade segment to the r = . 75R location. The next four 
terms are the first and second harmonics of cyclic pitch and 
(tanCj)^ is the coupled effect of the flapping angle on 
the pitch angle. 



(16) 0 = (9.is +0'- AIScosH* - Hissing - A2Scos2^> 

- B2Ssin2M» - t an ( ^ ^ ) cos 



The local Mach number is calculated as U/a where a is 
the speed of sound and U is given in equation 17. 

(17) 0 = (UP 1, + UT l + UR i )'k’ 

The angle of attack, can now be calculated. Figure 1 
illustrates that «< = e>- <f> . Initially, the program requires 

estimated values for the inflow ratio, collective pitch at 
the PSI equal zero position, harmonic cyclic inputs, blade 
twist and initial flapping angle and rate. The user can 
either input these values or accept the program's automatic 
default values of -.02, 5, -1.2, 7.53, 0, 0 and 0 

respectively. Using these assummed values, the initial 
angle of attack can be determined as fellows. 

UP = ^n.Rcos^ ♦ (r - E/R)^ + ( Vcos^ 5 ) cosYsin ^ 

reduces to 

UP = 7ML R 
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UT = (E/R + (r -E/R) cos^)-fL + { Vcosds,) sin ^ 



reduces to 



OT = (E/R + (r -E/R) ) -TL 



U R = (Vcos oi 5 ) ccs'tcos £ 
reduces to 

OR = Vcos 0^5 

In the program Vc os*( s is replaced by the term / £l 41R where/c 
is the advance ratio in the shaft axis system. Since the 
local pitch distribution and inflow ratio are estimated, the 
local angles of attack can be determined. The next section 
describes the method used for calculating flapping angle and 
rates at the N + 1 azimuthal position. 

C. METHOD OF SOLUTION FOR THE FLAPPING EQUATION 

In the preceding section, equation S was developed in 
order to calculate the time history of the flapping mcticn. 



differential equation for establishing the flapping angle as 
a function of PSI. The numberical solution is accomplished 
by use of a finite difference equation and a step-by-step 
procedure. An important characteristic of the solution is 
that it is periodic in nature. The function which 
represents a steady state flapping solution has the property 
£(H>) = ^ (v-Ht r) . Using this fact, a Fourier harmonic series 
can be written to describe the blade flapping motion. 

(18) ^ = AO - AlcosY - 3 Is i n^ - A2cos24> - E2sin2'+ > - 

A3cos3V - B3sin3^ .... 

3 v assuming that the first harmonic flapping is much larger 



* * 




The relation involves a complicated second order 
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than the 


other 


higher harmonics, the , 


series can 


be reduced 


to 










(19) 


§ 


= AO - AlcosSf - Blsin4» 






This 


equation can be differentiated with res 


pect to PSI 


to obtain 










(20) 




= AlsinM? - Blcos'j? 






(21) 


V 


= Alcos'V - B 1 sin*'? 

* 


* ^ 




Assuming that the values of p , (6 


and ft are known at 


some azimuthal 


position, the following 


equations 


must hold 


(22) 


Mg 


= AO - Alcosy - Blsinu) 






(23) 




= Alsin'V - Blcos^ 






(24) 


V. 


= AlcosV + Blsin^ 






Th ese 


equations can be solved for 


the N + 


1 azimuthal 


position 


by 


substituting Wa-m = Mfn * 


k'V into 


the above 


formulas . 


The 


flapping angle equation 


becomes 





(25) = AO - A1cos(%+^) - El sin ( 'V* +k$) 



By using the following two identities equation (25) can 



be rewritten 


as 


eg 


uati 


on (28) . 




(26) 


COS 


(^A + 


W) 


= cos^cosAH 1 


- s in^ sin 


(27) 


sin 


(^ + 


W) 


= sin'^cosi^'V 


+ cos'V A sinA.^ 


(28) 




= 


AO 


- A 


1 ( cos% COS^Vf 


-si nSw sinA.'V) 








- 


B 1 (s 


inH^cosA.y - c 


os^ sinA.vf) 


The 


terms 


c 


an 


be r 


earranged in 


to equation 29. 


pr ocedur 


e can 


b 


e u 


sed 


to develop eg 


uation 30 for ^ 


(29) 




= 


AO 


- c 


os ( A 1 cos 


+ 3lsinHV>) 




■ik 


+ 


si 


nW ( 


A1sinVr\ - Blc 


OS^rx ) 


(30) 




= 


cos IM( 


Alsin^vs - 3 1 c o s 4U ) 






+ 


si 


n &<*( 


AlcosV + 31s 


in'ftO 


Substitution 


in 


to 


equa 


tions 22, 23 


and 24 reduces 
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u 

IA 



















expressions to the flapping equations used in the program. 
The user can either enter the value for kV or accept the 
programs automatic default value of 15 degrees. 

* * *** 



(31) 




— 6^ Cos **• Sikj 












(3 2) 




- (3^ + Svvj CV ^ 




While 


this 


integration scheme is net one of the 


standard 


methods 


used , 


it is very useful in obtaining 


periodic 


solutions 


for differential equations similar to the 


one used 



here. Notice that for small values of the trigncme trie 
expression can be reduced to an ordinary Taylor series. 3v 
assuming sin&v equals A.'f and cos equals one equation 29 
reduces to 



(33) 



(jL, = 



By assuming sin equals LW and ccsLV equals 
twe terms of the cosine series 1 - equation 



reduced to 



(34) 
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D. FORCES, MOMENTS AND RADIAL INTEGRATION 

The forces acting on a rotor blade may be found by the 
summation of the elementry forces along the span at any 
azimuthal position. The forces considered by the program 
are the resulting aerodynamic forces only, and are initially 
summed in the shaft reference axis system. The program 
computes forces in three axis systems. They are the (1) 
shaft, (2) control and (3) relative wind axis systems. 

The program radially integrates differently for lift and 
drag calculations. The drag is calculated from the hinge 
offset, E/R , to the rotor tip. The lift is calculated from 
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the hinge offset to the next to last rotor blade segment. 
The last segment is considered a blade "Tip Loss Factor" 
segment. It is assumed that the tip trailing edge vortices 
cause no lift to be produced in this segment. The normal 
procedure is to define this segment as the last three 
percent of the rotor blade radius. 

The first of the maximum 15 blade segments is 
considered the spar or cut out segment. This is defined as 
the area between the hinge offset and the point where the 
airfoil actually begins. If no spar data are entered, it is 
asssumed that this first segment produces no lift and drag 
is obtained by using the CD verse Alpha Tables for the first 
blade airfoil data deck entered. If spar data are entered, 
the first segment and all other segments designated spar 
segments will have the lift and drag characteristics of the 
entered spar data. 

E. AZIMUTHAL INTEGRATION METHOD 



r oto 


r 


forces 


an 


d moments 


n in 


te 


grat ion 


a 


round the 


er 


to 


obtain 


JL 


he average 


solu 


ti 


cn of the 


f la p ping 


ep - 


by 


- step 


me 


thod, the 




are known 




only at a 


poin 


ts 


around 


th 


e azimuth. 


r al 


is 


a periodic function 


at io 


n 


can be 


sh 


own tc be 


s . 


Th 


is resul 




makes the 



azimuthal integration simple. The following method is used 
where N is the number of azimuthal positions used in the 
calculations and b is the number of blades. 

^n- K Ji K 

(35) JL. j \ \ V 

00 i - \ \® ' 
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F. MAJOR ITERATION 

Once the program has calculated a steady state flapping 
solution and has determined the resulting forces and 
moments, a question remains to be answered. Is this the 
desired solution and, if not, what must be done to obtain 
this solution? In order to solve the flapping equations 
certain known or assumed values were used, including inflow 
ratio ( 'A ) , collective pitch (9.75) and the cyclic pitch 
(A IS, BIS, A2S and B2S) . One method, first described in 
Ref. 3, is to iterate on the required lift and drag through 
modification of "A and 9.75, where 6.75 is the pitch angle 
at the 75 percent radius station at the PSI equal zero 
azimuthal position. The modified values are then reentered 
into the flapping routine and the calculation is repeated 
until it converges to within a specified tolerance on the 
required lift and drag. Drag is used here in the sense of 
negative rotor propulsive force. The program procedure is 
outlined below. 

The required rotor lift and propulsive force are 
expressed in terms of the magnitude and direction of the 
resultant force in the longitudinal plane. These are shown 
in Figure 6 where it can be seen that 

(3 6) a" = o( s + a' 



The required lift and propulsive force and their 
resultant are shown in Figure 7. In a similar way. 



(37) R^ = ( (F^f + (F^) Z ) ,h_ 



(3 8) L = R^cos (a ") 



D = Rj^ sin (a") 




•h. 



(39) R 



(40) a" = arctan (D^/L^) 
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The differences between the required and the computed R^ 
and a" values are defined as 



(4 1) \ 



(4 2) & a" = a£ - a" 

In order to correct 7\ and 9.75 to compensate for the 
difference between and the required 

values are expanded in a Taylor series with A and 9.75 as 
variables. The first order equations are: 



(43) 



= Rjl + ^ 



( 44 ) a.'L =. a! + 0 \ ** &©. n s- 

* ** 

Solving the equations for the iteration on and 9.75 

yields 



(4 5) = 



lal' 

^©•iS 



k<x!' 

^ ©OS 



Ha W' 

A ^.©n5 



IQd A a" 

©.T 'S' 



(“6) k© ^ V* 



^ A. a." 



v u 



^ W* _ 



X©.-»s 



How that R 4*.' and a n t and a" are known, the corrected 



values of 'X and 9.75 can be approximated by: 



(47) 



~/V , ~ -* r k.'A 



M 









- ©. 









In order to solve equations 45 and 46, the values of the 
partial derivatives in these equations must be found. The 
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Figure 6 



Resultant GRP Force Calculation Diagram 
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Figure 7 Required GRP Force Diagram 
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procedure used in based on the Wheatley-3ailey method and 
the formulas can be found on pages 186 and 207 of Ref. 4. 
Reference 5 outlines the derivation and a complete 
derivation was performed and verified in conjunction with 
this thesis. 



G. FLOW CHART 

The flow chart in Figure 8 is a block diagram showing 
the relationship between the various parts of the GRP 
program. 
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a solution. They are as fellows. 
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Normal Routine 



1 . 



In the normal solution to the problem the computer 
will vary the inflow ratio, LAMBDA, and the pitch angle, 
9.75, in order to produce the required lift and drag. The 
. cyclic inputs, A IS and BIS, are considered constants. 
Uniform inflow is assumed unless the user induces a 

non-uniform inflow by the use of variables 117 and 118, LAML 
and UVL. The program will calculate the required shaft axis 
angle and position of the control axis from the fixed value 
of BIS. In all of the options, the user has a choice of 
using either a program calculated first estimate of flapping 
angle, velocity and acceleration or an initial value of zero 
for all of the above. Either method usually requires the 
same approximate amount of computer time. The variable 
PCNV, item 97, determines the method tc be used. 



2. Des i re d Flappi ng Angl es 

This option allows the user tc specify the desired 
longitundinal and lateral flapping angles with respect to 
the rotor shaft axis. The program at each intermedient 
force iteration will vary AlS, BIS, inflow ratio and pitch 
angle. Variables 100 and 112, 113 and 114 control the use 
of this option. 

3 . Sho r t Iteratio n Sc h em e 



This option follows the 
major exception. The program wi 
through the flapping routine on ea 
program may or may not arrive at a 
flapping in the first few iterati 
the first few iterations in the 
rough estimates of the way the var 
and that an exact flapping solution 
if the user is only interested 



nor 


mal routine 


with 


one 


11 


make only 


one 


pass 


ch 


major iteration. 


The 


ste 


ady state sol 


utio 


n for 


cns 


It is assumed 


that 


nor 


mal routine 


are 


only 


iab 


les should b 


e ch 


anged 


is 


not actually 


reg 


uired 


i 


n transient 


£ la 


pping 
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behavior. In order to use this routine, the user must input 
a negative number for the variable XITLIM, item 73. The 
absolute value of XITLIM will still determine the maximum 
allowable number of times the program will enter the major 
iteration routine. 

4 . T ri mmed M o ment s 

This option will follow the Normal method described 
in paragraph one. The Normal method will only converge on 
lift and drag and will not consider the moments produced. 
If the variable TRIM, item 160, were assigned a non-zero 
value, the program would attempt to trim out the pitching 
and rolling moments about the rotor shaft. It is suggested 
that this be done in a two case run. The first case should 
be a normal run, with the desired printout. Then, for the 
second case, input the variable TRIM. This will do two 
things. First, it will provide a converged solution without 
consideration of moments. Secondly, it will allow the full 
number of iterations to be used to reduce the moments. Tne 
program does this by a short routine varying A IS and 31S. 
During this process, the whole flapping and force iterations 
must be repeated, but it will at least start by using a 
converged solution. The variable PCNV, item 97, should be 
non-zero to allow the flapping solution from the previous 
case to be used as a first estimate in this case. Setting 
the variable SKIPIN, item 91, equal to zero will allow a 
force and moment summation to be' outputed for each major 
iteration. This will allow the user to see exactly how the 
program is proceeding. 

5 • TO P Opt ion 

In this option, the program iterates upon the 
required lift but ignores any values inputed for required 
drag. This option can be used to simulate a wind-tunnel 
test. The user must input the shaft angle, item 111, and 
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the pitch 


a ngle , 


item 87 . 


These two 


inputs 


will be 


held 


co nstant . 


The program will 


iterate on 


the 


inflow ra 


tio. 


LAMBDA, in 


order 


to obtain 


a solution. 


The 


variable 


TOP, 


item 96, controls 


the use of 


this option 


• 







6 . AL OPT O ption 

This option, like TOP, is a wind tunnel option. It 
also iterates on required lift only. However, here the 
shaft angle, item 111, and the inflow ratio, item 88, are 
required inputs and are held as constants. The program will 
iterate on the pitch angle, item 87, in order to determine a 
solution. The variable ALOPT, item 110, controls the use of 
this option. 

B. GRP DATA DECK ORDER 

Data is entered into this program in the following 
or der . 

1. Airfoil Lift Coefficient Table 

2. Airfoil Drag Coefficient Table 
(Repeat steps one and two as necessary.) 

3. Case Input Data 

4. Harmonics of the Inflow Ratio ** 

5. Spar Lift Coefficient Table ** 

6. Spar Drag Coefficient Table ** 

** Optional Data 



C. AIRFOIL ELADE DECK REQUIREMENTS 



The GRP program requires that all 
be entered into the program in tabula 
to 15 different Mach numbers and five 
be entered. The program currently do 
values for the Coefficient of Moment 
segments of the retreating blade are 
region, angle of attack tables are 



CL and CD information 
r form. Tables for up 
different airfoils can 
es not use cr require 
, CM. Since certain 
in the reverse flow 
required to include 
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valued from -180 degrees to +180 degrees. If they are not 
included, an error message will be printed when the program 
can not locate a value for CL and CD at these large positive 
and negative angles of attack. If complete angle cf attack 
information is not available for a particular airfoil, the 
user can use the values provided in Section 2. It is 
realized that available data on airfoil behavior at large 
angles of attack are very limited, but so is the region on 
the rotor disc where the blade operates at these high 
angles. Since this occurs only immediately around and 
within the reverse flow region where dynamic pressure is 
low, little precision is lost in performance calculation by 
using one common representation for most airfoil behavior. 

As a minimum, two values for CL and CD at two different 
angles of attack and Mach numbers must be supplied. As a 
maximum, 15 different values of Mach number may be entered. 
Each Mach number may contain up to 48 different values of CL 
and 48 different values of CD and their associated angles of 
attack . 

The first Mach number must be equal to zero. This table 
can be an exact duplicate cf the lowest Mach number table 
the user has available. The highest Mach number table 
should be high enough in order to prevent the program from 
stopping because of a local Mach number higher than that in 
the table. A quick check can be obtained by adding together 
the rotor tip velocity and the forward flight speed. This 
combined velocity, divided by the local speed of sound, must 
be less than the maximum Mach number entered into the 
program. The program linearly interpolates between Mach 
numbers and angles of attack in order to determine the value 
of CL and CD. The subroutine 3LIN4 does the interpolation. 

Several options are available to help reduce the number 
of data points that must be entered. If the airfoil is 
symmetrical, the user only needs to enter values for 
positive angles of attack. The program will assign 
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the 



appropiate sign to the value of CL and CD according to the 
sign of the angle of attack. This is accomplished by case 
input variable 107, SYM. This option can also be used for a 
cambered airfoil where values of CL and CD at negative 
angles of attack are unknown. 

Values for large angles of attack need not be entered 
for each Mach number by making use of the program's input 
variables 156 and 157, or their automatic default values. 
Values for large positive and negative angles of attack need 
only be entered for the two lowest Mach number tables. The 
lowest Mach number table must be at a Mach number equal to 
zero. If an angle of attack is greater than variable 156 
(HIALFA) or lower than variable 157 (LOALFA) , or the 
programs default values of plus and minus 30 degrees, the 
Mach number is set equal to zero. This ensures that only 
the first two Mach numbers have to carry the whcle range of 
angles of attack from -180 to +180 degrees for a cambered 
airfoil or from zero to +180 degrees for a symmetrical 
airfoil. 

The format of the table input will now be described. 
The first data card contains the variable WELADE. This 
controls whether or not the user receives an echo printout 
of the Blade Airfoil Data being entered. If the value of 
W3LADE is equal to zero, the user will not receive an echo 
printout of the Blade Data. If W3LADE is a ncn-zero number, 
the user will receive the echo printout. The read format 
for WBLADE is F10.0. WBLADE is the only item on the first 
data card. 

The next card also contains only one piece of data, 
N3LADE. NBLADE is the number of different airfoil data sets 
to be used and appears on this card in an 12 format. This 
program will accept up to five different blade airfoil data 
sets. It will also accept one blade spar data set. If the 
rotor blade being analyzed were composed of three different 
types of airfoils, M3 LADE would equal three. However, if 
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the blade consisted of three sections, of which the first 
and third section were the same, N3LADE would egual two. It 
is explained later how the blade segments are assigned their 
respective airfoil type. This arrangement provides the user 
with the ability to vary the make-up of the blade while only 
having to enter into the program once a particular set of 
airfoil data. 

The above two variables, WBLADE and NBLADE, are only 
entered once for each complete computer run which uses the 
same set of airfoil data. The following information will be 
entered twice for each type of airfoil used. It will be 
entered first for CL and secondly for CD for each type 
airfoil. The overall format is summarized below. 

Card 1 W3LADZ 

Card 2 NELADE 

Card 3- CL's for airfoil number one 

CD's for airfoil number one 

CL's for airfcil number two 

CD's for airfcil number two 

Repeat as necessary. 

Card number three contains the variable NZ, which is the 
number of Mach numbers for which CL's will be entered for 
the first airfoil. This number must be right justified in 
12 format. The maximum number of Mach numbers for each CL 
and CD for one airfoil is 15. This is the only number 
entered on this card. 

Card number four begins the actual Mach number, CL 
versus angle of attack tables. The format in this paragraph 
must be repeated for each Mach number. This first card is 
divided into 11 fields (12, 10F7. 0) . The first field is a 
two-digit, right- justified integer in 12 format. It is 
egual to twice the number of data pairs for this Mach number 
plus two. This tells the computer how many numbers are 
reguired to be entered for this particular Mach number. A 
data pair consist of one angle of attack and its associated 
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CL or CD. The remaining ten fields are each seven columns 
long in floating point or F7.0 format. These fields begin 
in columns 3, 10, 17, 24, 31, 38, 45, 52, 59, and 66. On 
this first card, the field that starts in column three 
contains the number of data pairs at this Mach number. The 
field that begins in columns 10 is the actual Mach number. 
The remaining eight fields on this card are for the first 
four data pairs starting with the lowest angle of attack and 
increasing towards the highest angle of attack. If a 
symmetrical airfoil option is used, the lowest angle of 
attack is zero. If a non- symmetri cal airfoil is used, the 
lowest value for the first two Mach numbers should be -180 
degrees and for all the remaining Mach numbers the value of 
LOMACH entered or -30 degrees. 

All the remaining cards for this particular Mach number 
will contain the data pairs. These cards contain ten fields 
each, the first field consisiting of nine columns and the 
remaining nine fields consist of seven columns beginning in 
column ten and following the same format as the first card. 



The format 


is (F 9 . 0, 


9F7 


. 0) . This 


card is repeated as o 


ften 


as needed. 


Colu mns 


73- 


80 are not 


read and may be used 


for 


comments. 












This 


procedure 


is 


repeated 


until all the CL's 


are 


entered. 


Once completed 


the prog 


ram is ready to enter 


the 



values for CD. The whole procedure is repeated again, 
starting with the value for NZ representing the number of 
Mach numbers for which CD's will be entered. If more than 
one airfoil data deck is to be used, the above procedure 
will start over again by reading in the value of NZ for the 
number of Mach numbers to be entered for values of CL for 
the second airfoil. The number of data pairs for each Mach 
number must be between 2 and 48. The number of "twice the 
data pairs plus two" must be between 6 and 98. 

Prior to entering the airfoil data, the 
review the following input variables. 



user should 



and nonsymmet rical 



1. SYM (107) - Symmetrical 

airfoil data input control. 

2. HIALFA (156) - The highest angle of attack for 

which values of CL and CD will be found at all 
Mach numbers. 

3. LOAFLA (157) - The lowest angle of attack for 

which values of CL and CD will be found at all 

Mach numbers. 

4. SPAR (103) - Number of segments using spar 

airfoil data. 

5. TIPSWP (158) - Amount of tip sweep in degrees. 

6. TPS W ST (159) - Blade segment number at which 

the tip sweep begins. 

7. BSPL (120) - Input control variable for spar 
data . 

8. RB (I) (161-175) - Controls the blade segment 

airfoil data assignment. 



D. SPAR DATA DECK REQUIREMENTS 

The format for spar data are similar to that of the 
airfoil data with the exception that only one set of spar 
data can be entered into the program. Before the program 
will read spar data, input variable number 120, BSPL, must 
have a non-zero value. In addition, variable 103, SPAR, 
must indicate the number of blade segments which are using 
the spar data. The program automatically assumes that the 
first segment is a spar segment. This is further explained 
in the Case Input section. The spar data are the last tc be 
entered into the program. This is an optional input and is 
not required. If spar data are not inputed, the program 
will assume that the one automatic spar section creates no 
lift and has the drag characteristics cf the first airfoil 
section entered. 
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The format for inputing spar data are as follows, 
first card contains the variable WRSPAR in F10.0 form 
A non-zero value of WRSPAR causes an echo printout of s 
data to occur. The remaining spar data are handled 

exact same way as the airfoil section data, starting w 
the variable NZ. There is no input similar to that of 
airfoil section stating how many different spar data de 
are being entered since only one is allowed. As before 
minimum of two Mach numbers are required to be entered, 
blade segment printout indicates spar segments by the use 
a "0" for that segment. 

Input variables associated with SPAR data are 
follows. 



1 . 

2 . 



SPAR (103) - The number of blade segments using 
spar data. 

BSPL (120) - Input variable which controls the 
input and use of spar data. 



S. SAMPLE BLADE INPUT 

The next several pages illustrates the sample format 
a blade which has the following characteristics: (1) 

echo printout is not required, (2) there are two airf 
decks to be read in and (3) the first airfoil deck has n 
values of Mach numbers for which CL’s are to be entered. 



The 
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INPOT CASE LISTING 



GRP CASE INPUTS 



ITEM 

NO. 


DESCRIPTION 


PROGRAM 

VARIABLE 


DIMENSION 


REQUIR 


1 


Tip Speed 


OMEGAR 


EPS 


YES 


2 


Ra dius 


R 


FT 


YES 


3 


Speed cf Sound 


SPSD 


FPS 


YES 


4 


Air Density 


RHO 


SLG/CUFT 


YES 


5 


No. of Blades 


XNB 


- 


YES 


6 


Forward Speed 


VEL 


KTS 


89, 90 


7 


Offset Ratio of 
Flap Hinge (e/R) 


ER 


- 


YES 


8-22 


Delta X 


DX ( 1 5) 


- 


YES 


23-37 


Local Twist 


TW (15) 


DEG 


92 


38-52 


Local Mass Density 


XMASS (15) 


SLG/FT 


78,79 


53-67 


Local Chord 


C 


FT 


YES 


68 


Delta PSI 


DPS I 


DEG 


** 


69 


Flap Iteration Limit 


FTRL 


- 


** 


70 


Initial Beta 


BIN 


RAD 


** 


7 1 


Initial Beta * 


SPIN 


RAD/SEC 


** 


72 


Initial Beta ** 


BPPIN 


RAD/SEC**2** 


73 


Lift and Drag 
Iteration Limit 


XITLIM 


- 


* * 


74 


Required Lift 


RL 


LB 


YES 


75 


Required Drag 


RD 


LB 


95 


76 


Lift Tolerance 


XLTOL 


LB 




77 


Drag Tolerance 


XDTOL 


LB 


** 


78 


First Moment about 
Flap Hinge (M) B 


FMOM 


SLG-FT 


38-52 


79 


Second Moment about 


SMOM 


SLG-SCFT 


38-52 


80-82 


Shaft Orientation 


AG, BGL , CG 


DEG 


** 


83 


Pitch-Flap Coupling 
Angle (Delta 3) 


TD3L 


DEG 


- 


84 


Drag Increment 


DELD 


- 


- 


85 


Lat. Cvclic Pitch 


A 1 S 


DEG 


* * 
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86 


Long. Cyclic Pitch 


BIS 


DEG 


** 


87 


Collective Pitch 


T75 


DEG 


** 


88 


Inflow Ratio 


LAMBDA 


- 


** 


89 


Advance Ratio 


MUL 


- 


6, 


90 


VEL Control 


UIN 


- 


6, 


91 


Iteration Output 


SKIPIN 


- 


** 


92 


Linear Twist 


TWIST 


DEG 


23 


93 


No. of Blade 
Segments 


XNSEG 


- 


** 


94 


Climb Rate 


RCFPM 


FPM 


- 


95 


Flat Plate Area 


FPAREA 


FT 


75 


96 


Thrust Option 


TOP 


- 


- 


97 


Flapping Re-Use 
Indicator 


PCN V 


- 


** 


100 


MU Iteration ^ 

Tolerance 


ABIT 


— 


- 


101 


Beta Tolerance 


BTOL 


- 


☆ * 


102 


Beta* Tolerance 


3PTOL 


- 


* * 


103 


Spar Segments 


SPAR 


- 


- 


104-105 


Second Harmonic 
Control Inputs 


A2S, B2S 


DEG 


- 


1 06 


Solidity 


RSL 


- 


- 


1 07 


Symmetric 

Airfoil 


SYM 


- 


- 


108 


Spring Constant 


SFH 


FT-LE/RAD 


- 


109 


Damping Constant 


FDMP 


FT-LB/ 

RAD/SEC 


- 


1 10 


Lift Only Option 


ALOPT 


- 


- 


lu- 


Shaft Angle 


ALL 


DEG 


- 


ll 2- 1 1 3 


Desired A1 and 
B 1 Flapping 


RAIS 
RBI S 


DEG 


— 


1 14- 


A 1 , B1 Tolerance 


TOLAB 


DEG 


** 


1 15 


Tangent Delta 3 


TD3B 


- 


- 


1 16 


Phase Angle for 
Delta 3 


PRD 3 E 


DEG 


- 


1 17 


Induced Velocities 


LAM L 


- 


- 


1 18 


Induced Velocities 


UVL 


- 


- 


1 19 


Mot Used 








1 20 


Input Spar Data 


BSP L 


- 


- 
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1 21 


Azimuthal Printout 
Indicator 


PPSI 


DEG 


** 


1 22 


Minimum Lift 
Curve Slope 


ATEST 


- 


- 


123 


Iteration Gain 
Factor 


IGC 


- 


** 


1 24-125 


Not Used 








1 26 


Pre-Coning Angle 


PCR 


RAD 


- 


127-137 


Not Used 








1 38 


Hub Moment Inplane 
Aero Forces 


INPL 


- 


- 


1 39-1 41 


Aircraft Yaw, Roll 
and Pitch Angular 
Velocities 


PSIS, PHIS 
THFS 


RAD/SEC 




1 42 


CG Station 


FSCG 


INCHES 


- 


1 43 


Rotor Center Station 


FSMR 


INCHES 


- 


1 44 


CG Waterline 


WLCG 


INCHES 


- 


1 45 


Rotor Waterline 


WLMR 


INCHES 


- 


1 46 


Aircraft Lateral 
Velocity 


VELY 


KTS 


- 


1 47 


Spar Symmetry 


SYMSPR 


- 


- 


156 


Airfoil Tables 


HIALFA 


DEG 


* * 


157 


Airfoil Tables 


LOALFA 


DEG 


* * 


1 58 


Tip Sweep 


TIPS WP 


DEG 


** 


1 59 


Sweep Station 


TPSWST 


- 


** 


1 60 


Rotor Moments 


TRIM 


- 


- 


161-175 


Rotor Blade 
Airfoil Data 
Assignments 


RD 




* * 




** Program has automatic default 


value. 
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GRP INPOT DEFAULT VALOES 



ITEM 

NO 


PROGRAM VARIABLE 


DEFAULT VALUE 


68 


DPS I 


15.0 


69 


FTRL 


1 5.0 


73 


XITLIM 


15.0 


76 


XLTOL 


100. 


77 


XDTOL 


50.0 


81 


BGL 


90.0 


85 


A IS 


-1.2 * 


86 


BIS 


7.53 * 


87 


T75 


5 . 0 


88 


LAMBDA 


-.02 


91 


SKIPIN 


1 . 0 


93 


XNSEG 


15.0 


97 


PCN V 


1 .0 


101 


BTOL 


. 001 


102 


BPTOL 


. 001 


114 


T0LA3 


0.25 


121 


PPSI 


DPSI 


122 


ATEST 


5.0 ** 


123 


IGC 


1 . 0 


156 


HIALFA 


30. 


157 


LOALFA 


-30. 


159 


TPSWST 


16. 


161-175 


RB 


1 .0 




* 0. if ABIT, item 


100, is non-zero 




** -50 is TOP, item 


96, is non-zero 
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H. CASS INPOT DATA 



Many of the case inputs are self-explanatory by their 
name listing alone, however, many are not. This section 
will explain the input variables and case options available 
to the user. 

The program enters all the variables into a 200 element 
array called V (I) . Prior to entering case data, the program 
will automatically do two things. It will initialize the 
array V(I) to a value of zero. It will assign the default 
values listed in the previous section to those particular 
variables . 

The V (I) array is associated with the variable names by 
equivalent statements. The user needs only to enter values 
for the variables that are different from the default or 
initialized values. In a multiple case computer run, the 
user need only enter variables that are different from the 
preceeding case. If no new value is enter for a variable, 
the value from the previous case is carried over. 

1 . I tem 6 - Velocity - VBL 



one of two ways. The user will input either V2L, item 6, in 
knots or advance ratio MOY, item 39. The value of UIN, item 
90, determines which variable will be used. If (JIN is zero, 
VEL will be used. If (JIN is non-zero, M(JL will be used. If 
VEL is used, the program calculates the advance ratio by the 
following expression. 



The computer program will accept forward velocity in 




If M(JL is used , the program will calculate t he f 1 igh 
velocity by the following expression. 



. ^ -v*n 

M - * >> ^ -n. K \ 



2 . Item _l_z Hinge Offset - ER 

The offset ratio of the flap hinge, E/R, is the 
distance from the center of the rotor shaft tc the vertical 
flapping hinge, normalized by the rotor radius, item 2. 



3 • Item s 8-22 - Delta X - DX 





Delta 


X 


is 


the 


non-dimen: 


sionalized 


w idt h 


of 


each 


in dividual 


blade 


segment 


starting 


with segm 


ent 


num 


ber one. 


There 


may 


be 


up 


to 


15 


segments 


entered. 


The 


nu 


mber 


of 


wi dths 


ent 


er ed 


here 


must 


equal the 


value of 


Item 


93 


, XNSEG. 



XNSEG is the number of segments intc which the blade is 
divided. This can range from two to fifteen. It is 
recommended that a value of ten or more be used for XNSEG. 
If XNSEG were equal to 12, values of Delta X would be 
entered for items 8 to 19 and no values would be entered for 
items 20 to 22. The sum of ER plus the summation of the 
Delta X’s must equal one. ER is the non-dimensional width 
between the rotor shaft and the flapping hinge offset. Item 
number 8, which is the first segment width, represents the 
width between the flapping hinge offset and the point where 
the actual rotor blade airfoil begins. This area is known 
as the spar or cut out segment if no spar data were entered. 
If item 103, SPAR, is zero, the program will assume that 
this first section creates no lift and has the drag 
characteristics of the first inputed airfoil data section. 
Since this area experiences relatively low dynamic pressure, 
the calculations in this segment do net have an appreciable 
effect on the outcome of the program. The area between the 
shaft and the hinge offset, SR, produces neither lift nor 
drag in the program's calculations. The last segment 



is 



considered a tip loss factor segment. The lift is assumed 
to be zero, but drag is calculated in the normal manner. In 
previous runs, the width of this section has been set egual 
to three percent of the rotor radius, or Delta X egual to 

.03. 



4 . It em s 23 -37 - T wi st - TW 

The program has two option 
twist into the calculations. If th 
item 92, TWIST, can be used. If 
the twist can be entered in items 2 
entered for item 92, the program 
and ignore all values entered for 
the center of each segment can be e 
23. If the blade contains ten seg 
be entered and no values for items 
A word of caution is neces 
twist option, item 92. The twist 
at the 75 percent chord point, 
assuming that the twist starts at t 
linearly out to the rotor tip. I 
airfoil started at the 25 percent 
value of -9 degrees, a value of 1 
degrees would have to be entered 
twist distribution to be calculated 

5. I teras_33 z 52_ I _XMAS S 



s fcr entering geometric 
e blade has linear twist, 
the twist is non-linear, 
3-37, TW. If a number is 
will assume linear twist 
TW. The local twist at 
ntered starting with item 
ments, items 23-32 would 
33-37 would be entered, 
sary regarding the linear 
is considered to be zero 
The twist is calculated 
he rotor shaft and varies 
f the actual rotor blade 
radius point with a twist 
inear twist egual to -12 
in order for the correct 
by the program. 



The program provides two methods fcr entering the 
local mass density or the moment of inertia information. 
The individual mass density for each section can be entered 
in items 38-52, or the First and Second Moment, FMOM and 
SMOM, about the Flapping Hinge, can be entered in items 78 
and 79. If data are entered for item 79, the Second Moment 
about the Flapping Hinge, the program will use the 
information provided by items 78 and 79. If variable 79 is 
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equal to zero, the program will calculate the First and 
Second Moments for items 38-52 and ignore any value entered 
for items 78 and 79. 



6 . Items 53 -57 - Chord 



The program has two methods for entering local chord 
data. The local chord at each segment can be entered 
starting with segment number one. If the chord is a 
constant chord, the amount of data to be entered can be 



reduced 


to only 


items 


53 and 


54. 


If 


items 53 


and 54 are 


eq ual , 


the program assumes 


that 


f h. € 


chord 


is constant 


throughout the 


radius 


a nd 


will 


igno 


re any 


informa tion 


entered 


in items 


55-67. 


Therefore , 


for 


constant 


chord, only 


enter values for 


items 


53 and 


54. 


If t 


he rotor 


solidity is 



not inputed in item 106, the program will compute solidity 
as follows. 



CT 



\o C . tjT 
TT 



7 • I t em_68_ z _DPSI 

DPSI , or Delta PSI, is the incremental azimuthal 
value by which the program advances in its blade flapping 
and. force summation routine. This number must divide evenly 
into 360 degrees. A minimum value of five degrees is 
permitted. It has been found that for most cases decreeing 
the value of DPSI below 15 degrees does not improve the 
accuracy but does .increase the computational time. As an 
example, the rotor horsepower required for one particular 
run was 1096 RHP for DPSI equal to 15 degrees, 1095 RHP for 
DPSI equal to 10 degrees and 1097 RHE for DPSI equal to 5 
degrees. DPSI has a default value of 15 degrees. 
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FTRL 



8 . Item 69 - 

FTRL is the maximum limit on the number of times 
that the program will enter the flapping iteration routine 
in. search of a steady state flapping solution. The program 
has an automatic default value of 15 iterations. The 
program usually arrives at a steady state flapping solution 
within three to four iterations. The only method by which 
the user can actually determine the number of flapping 
iterations required would be to use one of the debug output 
options. However, these options will create a huge amount 
of output and the user is cautioned about their use. 

9 . It em s_ 70^7 2 BINz 3PPIN 

The program has three options on how to assign the 
values for the inital flapping angle, velocity, and 
acceleration at the PSI egual zero position. The user may 
(1) input the values, (2) have the program itself calculate 
initial values, or (3) accept the default values of zero. 
Option three is the most used option. It has been 
discovered that there is little or no difference in solution 
time between option two and three. Variable number 97, 
PC NV , controls which option is to be used. If PCNV is 
non-zero, the program will use either option one or three. 
The program initially assigns values of zero to BIN, BFIN, 
and BPPIN before the initial case data are entered. PCNV is 
assigned the default value of one. Therefore, with no 
values entered by the user for items 70-72 and 97, the 
program will start with an initial value for flapping angle, 
velocity and acceleration at the PSI equal zero position of 
zero . 

If the user sets PCNV equal to zero, the following 
formulas are used for determining the inital values. These 
formulas are derived from page 194 of Ref. 4. 
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§ = $ + m ^ 






\ -V 

it 



- M- 



c *i* B - * M 



v - ^ 



The term )$ is referred to as the Lock number where 

V — Cos 9 ^ 

^v> 



In a run where more than one case is executed at a 
time, a non-zero value for PCNV allows the previous case 
values for flapping angle, velocity and acceleration tc be 
used as the initial estimate for these variables in the next 
case. If PCNV is equal to zero, the inital values will be 
calculated by the above formulas. Most cases are run by 
accepting the default value of one for PCNV. 



10. Item 73 - XITLIM 

After the program calculates a steady state flapping 
solution for its estimated values of inflow ratic, pitch 
angle and shaft tilt angle, the program determines the 
forces and moments generated by this solution. If the 
forces and, optionally, moments do not meet the required 
amounts entered by the user, the program will calculate new 
values to re-enter the flapping routine. XITLIM determines 
the maximum number of times the program will compare the 
calculated values to the desired values. The program has a 
default value of 15 for XITLIM. A majority of the 
solutions require approximately five iterations to converge. 
Once the XITLIM limit is exceeded, the program will stop and 
printout an "Exceeded Limit" statement. The program 
automatically prints as output the number of Major 
Iterations it uses in calculating the solution. 



c 7 

v. / • 



The sign of XITLIM controls the Short Iteration Option. 
If XITLIM is a negative number, the program uses this 
option. A complete describtion can be found in the section 
entitled Main Iteration Options. 

1 1 . I te m 7 4 - Re qu ired L ift 

The required lift, RL, should equal the actual 
weight of the helicopter that the rotor system is 
supporting . 

1 2 . Item 7 5 - Requ i red Drag 

The required drag can be entered in one of two ways. 
It can be entered directly in item 75 or it can be 
calculated by the program by entering the aircraft's flat 
plate area in item 95, FPAR2A. The drag is equal to the 
flight path force that the rotor system must overcome to 
sustain level flight at a certain velocity. Items 75 and 92 
can be entered either as positive or negative numbers and 
the program will provide the correct forward flight solution 
by assigning the proper sign value internally. If a value 
is entered for FPAREA, the program will ignore any value 
assigned to RD. The drag for FPAR2A is calculated as 
follows. 



If item 92 is not entered, the user must supply the 
value for the required drag by using the above formula. 
There are certain options for which the program ignores or 
does not iterate on drag. An example of this would be when 
the program is used to estimate wind-tunnel test results. 
If these options are desired, see variables 96, TOP, and 
110, ALOPT. 
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1 3 . Items 76-77 - Tolerance s 

XLTOL and XDTOL are the lift and drag tolerances. 
The program will iterate until both the lift and drag are 
within the tolerances given by variables 76 and 77. The 
program has automatic default values cf plus and minus 100 
pounds for lift and 50 pounds for drag. 

14. Items 78 -7 9 - Moments 

Information regarding the use of FMOH and SMOM can 
be located in paragraph five on XMASS. 

1 5 . Items 80-82 - Sh aft Axi s 

The gravitational or weight vector must be oriented 
to the shaft axis of the rotor system. Figure 9 shows the 
positions of these angles. The shaft can be oriented in any 
desired direction. The program will automatically assign 
the proper values for normal helicopter flight. The default 
values for AG, 3GL , and CG are 0, 90, and 0 degrees, 

respectively. This orients the shaft in the vertical 

direction for normal flight. If any other orientation is 
desired, the user must enter the appropiate values for items 
80 to 82. 



16. Item 83 - Delta 3 

Delta three inputs are controlled by variables 83, 
115, and 116. These are TD3L, TD3B, and PHD3D, 
respectively. If the flapping hinge is connected in such a 
manner as to cause the blade to change pitch due to 
flapping, this is referred to as a Delta Three hinge. Item 
83 is the pitch-flap coupling angle; 115 is the tangent of 
the Delta Three Bar; and 116 is the phase angle for the 
Delta Three Bar. The program reduces the pitch angle at a 
particular azimuth and segment by the quantity TD3 where 
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Figure 9 Resolution of Gravitational ^orce 
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TD3 = TD3L + TD3B*sin (PSI + 
Normally, variables 83, 115, and 11 



PHD3D) 

6 are zero. 



1 7 . Item 84 - Drag I ncr eme nt 

The drag increment is a value of delta CD that is 
added to the value of CD obtained from the airfoil input 
data decks. This is added as a roughness factor that 
naturally occurs on blades that are used on production 
aircraft. A value of .002 is normally used. The value is 
not added to the drag calculations for spar data. 

1 8 . It ems 85 -86 - Cycl ic Pitc h 



by 


four 


variables. 


The 


o th 


first and second 


ha 


es 


A 1 S 


and BIS, 85 


an 


puts. 


The variables 


A 2 


second 


harmonic inputs 


1 


inputs and the 


B 



the 



The lateral and longitudinal cyclic pitch is 
riables. The program allows 
id second harmonics of eye; 

BIS, 85 and 86, control ■ 
variables A 2 S and 32S, 104 
rmonic inputs. The A's cor: 
the lateral inputs and the B's correspond 
longitudinal inputs. The program has automatic default 
vales for A1S, BIS, A2S and B2S of -1.2, 7.53, 0 and 0 

degrees, respectively. The user may enter different values 
if desired. Unless other options are indicated, the program 
will keep these cyclic values constant throughout the run 
and will vary inflow ratio and pitch angle to obtain a final 
solution. The program will automatically calculate the 

postion of the shaft axis by the momentum theory and will 
use the value of BIS to 
control axis. Figure 1 
There is a dif 
the program to seek s 
lateral flapping. Thi 
and 112-114. If thi 
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determ 
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ne 


the 
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the 
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angles. The program also has an 
rolling and pitching moments. Once a 
that meets the lift and drag toleranc 
and 3 1 S are varied to reduce the mom 
controlled by variable 160. Normally 
the user not inputing values for ite 
100, 112, 113, 114 and 160. 

1 9 . Item 8 7 - Pitch Angl e 

Variable 87 is the initial va 
pitch at the 75 percent radius stati 
azimuthal station. The program has a 
degrees. This value is used only to 
The program, under the normal run o 
value in the process of iterating for 
There is an option where the pitch ang 
a wind-tunnel test. This is the TOP o 

2 0 . Ite m 88 - Infl ow Rat io 

The variable LAMBDA controls 
an uniform inflow. Since the equatio 
done in a gyrocopter mode, inflow is n 
down through the rotor. This is the 
mode. The program has a default valu 
The program will iterate on LAMBDA in 
rountine. The ALOPT, variable 110, o 
co nstant . 

2 1 . Ite ms 89-90 - MOL 

Information regarding the use 
found in paragraph one on Flight Veloc 

22. Item 91 - SKIPIN 

Information regarding the use 
in the section of Case Optional Output 



option to remove hub 
solution is obtained 



es. 


the 


val 


ues of 


A1S 


ents 


• 


This 


opt io 


n is 


r r 


uns 


ar e 


made 


with 


ms 


% 

in 

CO 


86, 


104, 


105, 



lue o 


f the 


collec 
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23. Item 92 



Linear Twist 



Information regarding t 

paragraph four on local twist. 

24. Item 93 - XNSEG 

Information regarding X 

paragraph three on Delta X. 

25. Item 94 - Rate of C l imb 

The program can be made 
solution for any given rate of cl 
descent rate must be entered in 
with positive values for climbs 
descents. The program assumes a 
across the entire rotor surface e 
or descent. This value is a 
correction into the calculations 
and angle of attack. 

26. Item 9 5 - FP AREA 

Information regarding the 
area, can be found in paragraph tw 

27. Item 96 - TO P 

The TOP option is one of 
If TOP is a non-zero number, the 
the required lift of variable 
required drag of variable 75. Ite 
at the 75 percent radius at PSI 
constant. Item 88, the inflow rat 
major iteration routine. A non 
result in a value of -50 for item 
minimum accepable value for the li 



wist can be found in 
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to calculate a complete 
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96 or 110 is executed. If non-zer 
ALOPT are entered, the program w 
Shaft angle, item 111, must te inpu 

28. Item 97 - PCNV 

Information concerning the 
Solution Re-Use Indicator, can be f 
Initial Flapping Conditions. 

29. Ite m 9 8 - PRINT 

Information concerning t 
program's main output indicator, ca 
entitled Case Optional Output Indie 

30. Item 99 - XE ND 

Item 99 is the End of Case 
last data variable that will be en 
XEND is a negative number, the pr 
determines a solution for that pa 
since an infinite number of cases 
computer run, XEND also tells the 
cases to go. If XEND is equal t 
assume that the next case will 
airfoil data. If XEND is any oth 
the program assumes that the next 
airfoil data and will enter only c 
the options which normally follow t 
time that the variable XEND is 
careful to follow the format for N 
is the number of inputs per data 
negative number for this XEND car 
sign on NNUM which actually key 
reading data cards for a particular 



o values for both TOP and 
ill do the TOP option, 
ted by the user. 

use of PCNV, the Flapping 
ound in paragraph nine on 



he use of PRINT, the 
n be found in the section 
ators . 
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31. Item 100 



ABIT 



Information regarding the use of A3IT can be found 
in paragraph 40. 

32. Item s 1 01-10 2 - B TOL 

In the blade flapping routine, the program searches 
for a steady state flapping solution for the given 
conditions of inflow ratio, pitch angle, and cyclic input. 
The program compares values of flapping angle and velocity 
at the PSI equal zero azimuthal position on each revolution. 
If at this position, the difference between the n-th and the 
(n - 1)th revolution values for flapping angle and velocity 
is less than BTOL and BFTOL, respectively, the program 
assumes that it has determined a steady state solution. 
BTOL and 3PT0L have default values of 0.000001 radians and 
radians per second, respectively. If other values are 
desired, the user may enter those values for items 101 and 
102 . 



33. Item_1_ 03_^_SP A H 

The number of blade segmen 
inputed in item 103. If nc spa 
value for SPAR should be enter 
program assumes that the first se 
the flapping hinge and the point 
begins on the rotor blade. This a 
the "cut out" segment. In this 
that this cut out area produces 
information from the first ai 
calculations. If a non-zero num 
spar airfoil data must be availa 
120, BSPL, controls the spar input 
to be inputed, BSPL should be assi 
a multiple case run, where spar d 



ts using spar data can be 
r data are available, no 
ed. For this case, the 
gment is the area between 
where the actual airfoil 
rea is also referred to as 
case, the program assumes 
zero lift and uses CD 
rfoil section for drag 
ber is entered for SPAR, 
ble. Case input variable 
option. If spar data are 
gned a non-zero value. In 
ata are initially entered. 
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the variable BSPL is set equal to zero as soon as the spar 
data are entered. If the program did not do this, the user 
would have to enter a zero for BSPL for the next case if no 
new spar data are to be entered. If in a multiple case run, 
the user decides to enter new spar data, the variable BSPL 
must be assigned a non-zero value for that particular case. 



34. Items 104-105 - A2 S B2S 



Information 
control inputs can 
Longitudinal Cyclic 



regarding the use of the 
be found in paragraph 18 
Inputs. 



second harmonic 
on Lateral and 



35. Item 106 - Solid ity 



Information regarding the use 
solidity, can be found in paragraph six cn 



of RSL, rotor 
local chord. 



36. Ite m 107 - SYM 



SYM is the non-sym 
If the user assigns a non- 
will assume that all blade 
The user must enter values 
range of angles of attack f 
value of SYM is zero, only 
degrees need to be entered 
variable 147, SYMSPR. SY 
exactly in the same manne 
data. 



metrical a 
zero value 
airfoil dat 
for CL an 
rom -180 to 
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37. Items 108-10 9 - SFH FDMP 



Values 

constant, FDMP, 
known. These 
hingeless rotor 



for the spring constant, SFH, and damping 
about the flapping hinge can be entered if 
variables can be entered to simulate a 
system or a system with flapping springs. 
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38. Item 110 



ALOPT 



This is one of the wind-tunnel options. If ALOPT is 
a non-zero input, the program will iterate to obtain the 
lift required of variable 74 but will ignore the required 
drag, variable 75. In this option, variable 87, collective 
pitch, will be varied, but not variable 88, inflow ratio, in 
the program calculations. This is the opposite of the TOP, 
varialbe 96, option. 

If a lift curve slope less than ATEST, item 122, is 
calculated while using this option, the program will stop 
and produce the following message. "Stall Criterion has 
been violated — will go to next case, if any." Item 111, 
the shaft angle, must be inputed. If non-zero values are 
entered for both TOP and ALOPT, the program will do the TOP 
Option . 

39. Item 111 - Sh af t A n gle 

Item 111, the shaft angle, must be inputed whenever 
the TOP or ALOPT options are used. 

40. Ite ms 1 12- 114 - RAIS 

If variable 100, ABIT, is non-zero, the program will 
iterate the blade flapping solution in an attempt to obtain 
the desired lateral and longitudinal flapping angles 
indicated by variables RAIS and 3B1S, respectively. The 
program will iterate to the accuracy indicated by TOLAB, 
item 114. Item 112 is RAIS and item 113 is REIS. 

4 1 . Items 115-_1_1 6_-_Del ta_3 

Items 115 and 116 are the tangent and phase angle of 
a Delta 3 Bar. Information regarding the use of TD3B and 
PHD3B can be found in paragraph 16 on the Pitch-Flap 
Coupling Angle. 
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42. Item s 11 7-118 - LAML 

The program allows the user tc induce a velocity of 
any form onto the rotor system. This is done in a harmonic 
series of the form (AO + A 1*cos (PSI) + 31*sin(PSI) + 
A2 *cos (2*PSI) + B2*sin (2*PSI) + ....) for each segment that 
the blade is divided into. If variable 117, LAML, is a 
non-zero number, the program will enter the harmonics of the 
induced velocities. During the Read routine, LAML, will be 
assigned a value of zero. Therefore, for each case where 
different values for the harmonics are desired, LAML will 
have to be set to a non-zero number. Variable 118, UVL, 
controls the use of the induced velocities. If UVL, is 
zero, the induced velocities will all be set equal to zero. 

A short example will now be given on the use of the 
control variables in a multiple case run. Assume that no 
harmonic induced velocites are desired for the first case. 
The user would make no inputs for LAML and UVL. For the 
second case, assume that harmonic induced velocities are 
desired. LAML and UVL would be set equal to a non-zero 
number. The harmonic variables would follow the case input 
data for this particular case. For the third case, it is 
desired that the same induced velocites be used. The user 
would not enter any values for LAML and UVL since (1) LAML 
has been automatically set to zero and hence nc new harmonic 
data will be entered and (2) UVL is still equal to a 
non-zero number. It is desired for case four to use no 
induced velocities. The user now would enter the value of 
zero for UVL and program will zero out all the induced 
velocity variables. 

This paragraph will describe hew to use the variable 
induced velocity option. The inputs A1, B1, A2, 92, .... 
are numbers in units of feet per second. If the direction 
of the velocity is down, a negative sign is associated with 
the A ' s and B's. A negative sign will decrease the angle of 
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attack for an element originally at a positive angle of 
attack. It will increase the amount of negative angle of 
attack for a blade element originally at a negative angle of 
attack. The values of A’s and 3's do not effect the 
uniform inflow ratio, LAMBDA, that the program iterates 
upon. The format for using the variable inflow velocity is 
as follows. 

1. A value for NHARM is entered in 13 format. 

NHARM is the maximum degree cf the harmonics 
that is to be used. If the maximum degree used 
is A3*cos (3*PSI) , then NHARM is three. 

The next three paragraphs are repeated for each 
blade segment. 

2. The value for AO is entered in E15.6 format. 

The I3M 360-67 at the NPS will accept F15.0 
format. 



3. 


The values for A1, A2, 


A3, . . 


. , up to NHARM are 




entered on the next 


data 


cards in 


format 




5S14 .6. 








4. 


The values for 31, 32, 


B3 , . . 


. , up to NH 


ARM are 




entered on the next 


data 


cards in 


format 




5E14 .6. 








5. 


Three cards are require 


:d for 


each blade 


segment 



whether or not the values are equal to zero. 

If the user has 15 blade segments, a minimum of 
45 data cards are required. 

43. It «S_120_^_BSPL 

Information concerning the use of 3SPL can be found 
in paragraph 33 entitled Number of Spar Airfoil 
Segments. 



Data 



44. Item_222_^_PPSI 



P P SI represents a delta PSI for printout purposes. 
PPSI has a default value equal to DPSI, variable 68. PPSI 
shall never be a smaller increment than the incremental DPSI 
used to calculate the solution. 

45. Item 122 - AT EST 

ATEST is the minimum acceptable value of the lift 
curve slope when option TOP or ALOPT is used. If no value 
is assigned, ATEST has a default value of 5 (1/rad) . Nher. 
the TOP option is used it has a automatic value of -50 
(1/rad). The lift curve refers to the increase in rotor 
lift with the tilting back of the tip path plane. A check 
on it for a minimum is for convergence purpose only. 

46. Item 123 - I GC 

IGC is the Iteration Gain Control factor. If the 
major iteration fails to converge, choosing a fractional 
value for IGC can greatly speed convergence. This may be 
especially helpful when parts of the rotor are in stall. 
The amounts by which the convergence algorithm changes the 
independent variables is multiplied by IGC. Setting item 
91, SKIPIN, equal to zero may help the user decide if this 
IGC option might be useful. 

47. Item_J126_ z _?CR 

The pre-coning angle in radians may be entered here. 

48. Item_J38_-_INPL 

Input 1.0 for IN PL in order to remove hub moment 
inplane aerodynamic forces from the calculation of 
aerodynamic pitch and roll moments about the shaft axis. 
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49. Items 139-141 - PSIS 

Variables PSIS, PHIS and THFS represent aircraft 
angular velocites. The aircraft can be given any angular 
velocity in yaw (PSIS) , pitch (PHIS) , and roll (THFS) in 
radians per second by the use of these variables. They 
effect the calculations of DP, UT and DR. Paragraph 50 
contains additional information. 

50. Items 14 2- 145 - CG 

If the angular velocities, items 139 - 141 are 

entered, they are used in the calculations of DP, UT and UR. 
If no values for items 142 - 145 are entered, the program 

assumes that the rotor system rotates due to the angular 
velocities about the center of the shaft. If values are 
entered for items 142 - 145, the calculations will assume 

that the entire rotor shaft is rotating about the center of 
gravity. Variables 142 - 145 are FSCG, FSMR, WLCG and WLMR . 
FSCG is the longitudinal CG position while FSMR is the 
longitudinal position of the main rotor. WLCG is the CG 
waterline station while WLMR is the main rotor's waterline, 
or vertical position. All of these values must be entered 
in units of inches. 

5 1 • I te m_J.46_^_V FL Y 

The aircraft's lateral speed in knots can be entered 

here. 

52. Item 147 - SYMSPR 

If the spar data are symmetrical, do not enter a 
value for SYMSPR. The program will only use values for CL 
and CD between zero to +180 degrees. If the spar data are 
no n- sv mm et r ical, enter any non-zero value for SYMSPR. The 
program will require values from -180 to +180 degrees. 
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53. Items 156-157 



HI ALFA 



In order for the program to p 
region in and around the reverse flew 
CL and CD are required to be known at 
attacks approaching 180 degrees. Howe 
the user from having to enter a who 
attack for all Mach numbers, the progr 
for angles above HI ALFA, 156, and b 

Mach number is set equal to zero for 
The user is only required to enter 
first two Mach number tables. HI A 
default values of +30 and -30 degress, 

54. Itens 1 58 -15 9 - Ti p S wee p 

The GRP pregram will properly 
sweep angle, UT, (JR, and pitch angle o 
TIPSWP is the amount of sweep of 
degrees. TPSWST is the blade segmen 
sweep begins. The program assumes 
outboard segments starting with TPSWST 
of degrees indicated by TIPSWP. DT a 
these segments as shown is Figure 11. 

55. Ite m__160_-_TRI M 

If TRIM is a non-zero valu 
attempt to adjust A1S and BIS in order 
and pitching moments to less than plu 
It will first obtain a solution wh 
required lift and drag, then it will 
reduce the moments. 

56. I tems_26_1-J_7 5_2_R B 

Variables 161 - 175 assign th 

to blade segments one through fiftee 
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adjust A1S and BIS to 



e proper airfoil data 
n. The program will 







Figure 11 Tip Sweep Calculation Diagram 
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accept up to five airfoil data decks. The RB array is 
initialized with the value of one for all 15 segments. If 
only one blade airfoil deck is used, there is no need to 
enter any values for RB. Only enter values for the segments 
which will use airfoil data sets two through five. If no 
spar data are entered, SPAR equals zero, segment number one 
is considered the cut out segment and is assigned the value 
of RB ( 1) equals one for calculation purposes and the value 
of RB ( 1) equals zero for printout indentif ication . In this 
case, segment one produces only drag tut no lift. If SPAR 
is greater than zero, RB(I) ,1=1, SEAR, will be assigned 
the value of zero for printout in dentif ication. These 
segments will use the spar data entered and will calculate 
both lift and drag. 



I. CASE INPUT FORMAT 



Below is located the format that is used to inpu- case 
data to the program. All input data cards are cf the format 

(12, 14, 5F12.0). The input cards contain seven fields 
which are called NNUM, NLOC, C(1), C ( 2) , C (3) , C(4) and 
C(5). An example of a data input card is as follow. 

5 1 700. 26.8 1148.6 .002246 4.0 

1. NNUM 



of 



NLOC 



NLOC 





is the 


numb 


er 


of i 


n 


put 


s on this 




card. 
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C(N) are 
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in the 
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Rotor Tip Speed. 
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3. 



4. 



C (1) 

to t 

Each 

in c 

31 - 

for 

for 

Radi 

of B 

Omis 

wit h 

grea 

Omis 

valu 

indi 

to 

200 , 
prop 
with 
err o 
fi el 



through C (5) are the val 



he in 


put items NLOC thro 


value 
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contain a 
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7 - 


18 for C (1) 


, 1 


42 for C ( 3) 


, 43 - 54 


for 


C (5) . 


In th 
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r D 
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J. CASE OPTIONAL OUTPUT INDICATORS 



The program has two variables 
printout, variable 91, SKIPIN, and 
will automatically produce a printo 
for each case and a one-page 
conditions and iteration limitatio 
the program. It will also produce 
resulting forces, moments and cal 
for the final converged solution, 
an echo printout of the airfoil a 
this echo printout is desired, 
correct printout indicators descri 



that control the ou 
98, PRINT. The pro 
ut of the case input 
summary of the in 
ns the user placed 
a one-page summary of 
culated rotor hcrsep 
The user can also rec 
r.d spar data decks, 
the user will find 
bed in sections C a 



tput 
gram 
data 
ital 
upon 
the 
ewer 
eive 
If 
the 
nd D 
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entitled Blade and Spar Data Deck Requirements. 

The main optional printout variable is variable 98, 
PRINT. PRINT can be a number from 1 to 1,111,111 depending 
upon the option desired. If no value is enter for PRINT, 
the user will receive the printout described above. Below, 
is listed the PRINT Options. If, fcr example, PRINT is 
assigned a value of 111, the user will receive printout 
options 1, 10 and 100. 



OPTIONAL OUTPUT INDICATORS 



Option 1 



Option 10 

Option 100 
Option 1000 
Option 10000 



Angle of attack, Mach number, 

section lift and drag 

coefficients, inflow angle, lift, 
and sweep angle at each azimuthal 
position for each radial blade 
segment. Only for converged 

flapping solution. 

Converged flapping angle, rate, 
and acceleration at each 
azimuthal position. 

Converged integrated forces on 
blade at each azimuthal station. 
Harmonic analysis of blade forces 
for converged case. 

Harmonic analys is of air loads 

for converged case. 



DEBUGGING OR TRANSIENT OPTIONS 
Option 100000 Transient flapping angle, rate, 

and acceleration at each 
azimuthal station. 

Option 1 C00000 Option 100000 plus blade 

velocities, angles, Mach number, 
section coefficients, and lift 
for each blade segment at each 
azimuthal station. 
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The user is cautioned that the debugging options 
can give a huge amount of output data. 

The second printout option is the variable SKIPIN, 
number 91. This variable controls the summary force, moment 
and horsepower output discussed in the first paragraph. If 
SKIPIH is greater than zero, this summary will be printed 
only for the final converged solution. If SKIPIN is equal 
to zero or a negative number, this summary will be printed 
for each loop through the major iteration (force summation) 
routine. SKIPIN has a default value of one. If the program 
is not converging to a solution, the user can see 
immediately, with very little extra printout, exactly what 
intermediate solutions the program is producing by setting 
SKIPIN equal to zero. This may help the user in deciding 
whether or not to use the Iteration Gain Factor, IGC, 
variable 123. 

K. IBM 360 EXECUTION CONTROL CARDS 

This section illustrates the control cards required to 
execute the GRP program using the IBM 360 at the Naval 
Postgraduate School. The program may be run under OS or 
CP/CMS. There are two ways of running the program under OS. 
The first is to run the entire program and data through the 
computer at the same time. The second way is to store the 
main program on a disk as a library program and enter only 
the data through the card reader for each desired case. The 
second method has the two advantages of (1) not requiring 
the user to enter the entire 1100 card main program through 
the card reader for each run and (2) the amount of CPU time 
required can be reduced since the main program does not have 
to be recompiled for each run. Ihe program requires 
approximately one minute and forty seconds of CPU time to 
compile. The normal run time for each case is approximately 
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20 CPU seconds. This can vary with the amount of printout 
data requested. If the program is compiled on the CP/CMS, 
the user will have to request 344K bytes core size on the 
login message. The standard 256K bytes core size is not 
large enough for compiling. However, once the program has 
been compiled, it can be executed within the 256K normally 
available. 

The following cards are required to execute the entire 
program through the card reader at one time. 

Standard Job Card 

// EXEC FORTCLG, REGION. FORT=150K 
// REGION. GO=180K 
//FORT . S YSIN DD * 

Main GRP Program 
/* 

//GO. S YSIN DD * 

Case Input Data 
/* 

The following two programs are used to reserve space and 
load the program onto a disk. 

Standard Job Card 
// EXEC PGM=IEFBR 1 4 

//LOAD DD DSN = S1395.HELO, UNI T= 3330, VCL=SER=DISK01 
// DISP= (NEW, KEEP) , LAB EL= RETPD= 1 50 ,S PACE= (CYL, (1,1,1)) 

/* 

Standard Job Card 

// EXEC FORTCL, REGION. FORT=180K 

//FORT . S YSIN DD * 

Main GRP Program 
/* 

// LINK. SYSLMOD DD UNIT=3330 , VO L=S ER=DI SK 0 1 , 

// DSN=S1395. HELO (GRP) ,DISP=SHR 

The S1395 used above and below must be change to S for 
student or F for faculty with the appropiate user number 
instead of 1395. The following cards must be used to 
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execute the program once stored on a disk. 

Standard Job Card 

//GO EXEC PGM=GRP, REGION=180K 

//STEPLIB DD 0NIT=3330,VOL=SER=DISK01,DISP=SHR, 

// DSN=S1395.HEL0 

//FT06F001 DD SYSO 0T=A , DCB= (RECFM=FBA,LRECL=133, BLKSIZE=3325 

//FT05F001 DD * 

Case Input Data 

/* 

J. Sample Program Output 

This section describes in detail the output available 
from the GRP program. In addition, a sample computer output 
is included with each describtion. The program will print 
up to ten different tables. Seven of these tables are 
optional and are not automatically printed. The ten tables 
are as follows. 

1. Echo Printout of Rotor Blade and Spar Data 

2. Case Input Data Card Listings 

3. Summary of Input Data 

4. Debugging or Transient Information (Options 
100,000 and 1 ,000 , 000) 

5. Summary of Forces, Moments and Horsepower 

6. Converged Flapping Solution (Option 10) 

7. Converged Integrated Forces (Option 100) 

8. Harmonic Analysis of Z Force (Option 1000) 

9. Converged Blade Analysis (option 1) 

10. Harmonic Analysis of Air Loads (Option 10,000) 

The input variable PRINT, item 98, controls the output 
of items four and six through ten. Items two, three and 
five are always outputed. Item one is controlled by the 
first data card on the rotor blade and spar section input 
decks . 
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1 . Echo Ro tor B lade Prin tout 

The next page contains a partial sample Echo 
Printout of Rotor Blade Input Data. This printout 
illustrates that (1) the printout was requested (WBLAD2 = 
10.), (2) there are two airfoil decks to be entered, (3) 
there are nine Mach numbers for which CL’s are to be read in 
and (4) the remaining portion of the printout is the values 
of Mach numbers, angles of attack and lift coefficients for 
the first airfoil deck. 
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2 • In pu t Data Car d Lis ti ng 



The next page contains a sample Case Input Data Card 
Listing. This is one of the automatic printouts. It is an 
echo printout of the input data cards 
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3 . Sum m ary of Inpu t Data 

The next page contains the automatic Summary of 
Input Data printout. The following information can be seen 
on this sample output. 

a. The blade was divided into 15 segments starting 
from the hinge offset and proceeding outward. 

b. No values were entered for the local mass 

desnsity, input variables 38 - 52. Instead, values for the 

First Moment, M (M) = 85, and the Second Moment, MOM 

INERTIA = 1450, about the Flapping Hinge, input items 78 and 
79, were entered. 

c. The X row indicates the calculated centers of 
each segment expressed in a percentage of the distance out 
the rotor blade. 

d. The BLADE DECK row indicates that the first 
blade segment was considered a Spar segment. Segments two 
through five and twelve through fifteen belong to airfoil 
data deck number one, while segments six through eleven 
belong to airfoil data deck number two. 

e. The rest of the information, with one exception, 
is a summary of the case input data. The exception is the 
term THRUST FACTOR. This is the value used to 
no ndimensionalize all the calculated forces in the program. 
The THRUST FACTOR equals TT C -0- fO . Moments are 
nondimensionalized by the THRUST FACTOR times the radius. 

f. The program checks to see if all the blade 
segments, delta X's, plus the distance from the shaft to the 
hinge offset, E/R, add up to one. If, on the printout, 
SUM(DX) + E/R does not equal one, the user has made a 
mistake somewhere with the delta X's or in the E/R number. 
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Deb uggin g or Trans ien t In fo rm a tio n 

The next page illustrates the Debugging or Transient 
Printout. This is PRINT Option 1,000,000. The information 
available includes flapping angles, rates and accelerations 
at each azimuthal position with a radial position display of 
OP, UT, U, PHI, Angle-of- Attack, Mach Number, CL, CD and 
Lift Per Inch produced. 

This option is generally outputed only when the user 
is experiencing unknown difficulties with the GRP program. 
The program will output all of the above information for 
every revolution and iteration until a converged solution is 
obtained or the program runs out of allowable computer time. 
If desired, the variable SKIPIN, item 91, will provide a 
printout of the Forces, Moments and Horsepower Summary after 
each major iteration. 
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VI. 



Fo rces S umma ry 



5. 



The next printout illustrates the Summary of Forces, 
Moments and Horsepower. This page is automatically outputed 
for the converged solution. A printout of this summary can 
be obtained for each loop through the major iteration 
routine by the use of the input variable SKIPIN, item 91. 
The following information can be observed from this sample 
pr intout. 

a. The cyclic lateral and longitudinal inputs, A1S, 
BIS, A2S and B2S are printed at the top of the page. This 
example indicates that there were no second harmonic inputs 
for A2S and B2S , which is normal. 

b. THETA .75 is the pitch angle of the blade at the 
75 percent radius station at the PSI equals zero azimuthal 
position. In most options the program will iterate upon the 
values for THETA .75 in its search for a converged solution. 

c. LAMBDA refers to the converged value for the 
uniform inflow ratio. If a rate of climb or descent was 
used in the case, that rate divided by the tip speed would 
have to be subtracted from or added to this value of LAMEDA, 
re spectively . 



d. MU(X)S and MU(Y)S are the advance ratios in the 
longitudinal and lateral directions with respect to the 
shaft axis. 

e. CT, CQ, CH, CL and CD are the calculated overall 
coefficients of Thrust, Torque, H Force, Lift and Drag. All 
of these items are nondimensionali zed by the value of the 
Thrust Factor. 

f. Lift and Drag forces are calculated with respect 
to the relative wind axis. Thrust and H forces are 
calculated with respect to the control axis. Z forces are 
calcuated with respect to the shaft axis. X and Y forces 
are calculated perpendicular to the shaft axis. 

g. The Equivalent Drag is the total drag force 
created by the fuselage, flat plate area times dynamic 
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pressure, plus the profile drag create 
system . 

h. The Equivalent ?. A., or 
area is obtained by dividing the Eg 
dynamic pressure. 

i. Alpha (S) is the shaft a 

Alpha (c) is the control axis orientat 
momentum theory to calculate Alp 
determined from the following rela 
Alpha (C) + BIS. 

j. LAT. DIS. and LONG. DIS. 
offset as a percentage of the rotor 
shaft to create the program's rolling 

K. PM and RM are the calcula 
and roll moments about the shaft axis, 
and Roll Moments are calculated 
aerodynamic, inertia and elastic hub r 



d by the turning rotor 

Equivalent Flat Plate 
uivalent Drag by the 

xis orientation while 
ion. The program uses 
ha (S) . Alpha (C) is 
tionship. Alpha (S) = 

are the lift vector 
radius from the rotor 
and pitching moments, 
ted aerodynamic pitch 
The SHEARS Hub Pitch 
from summing the 
estraint moments. 
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6 • F la p pin g So lutio n 

The next page contains a sample printout of the 
converged flapping solution. The first item tc appear are 
the flapping values at the PSI equal zero azimuthal position 
for each revolution prior to the converged revolution on the 
final iteration. The values for flapping angle and its 
first two derivatives are expressed in radians. 

The second part of this printout is the actual 
flapping values for the converged revoltion. The difference 
between the flapping angles and rates between the PSI equal 
zero and 360 degree position must be less than the 
tolerances entered in variables 101 and 102. The numbers 
here are also in radians. The last item is a Fourier 
coefficient series for the flapping angle and it is 
calculated in degrees. All calculations are done in 
respects to the shaft axis. This is PEINT Option 10. 
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7 • Force In tegratio n 

The next table is 
This is PRINT Option 100. 
blade only at a particula 
CQD are the Coefficients 
Torque due to Drag. CQ 
CX, CY and CZ are all re 
reference system. CMHS 
Moment due to aerodynamic 
moments about the Shaft A 
Rolling Moment due to th 
the (B) character is th 
example is four. SIG is 
MAX B *CQD/S IGM A is a 
calculated by determini 
between the PSI equal 1 
that has a value of CQD g 
PSI equal 180 degrees azi 



a sample Force Integration Output. 
The forces are a printout for one 
r azimuthal position. CQ, CQL and 
of Torque, Torque due to Lift and 
is calculated from CQ = CQD - CQL. 
lated to forces in the shaft axis 
is the Coefficient of Pitching 
, inertia and hub elastic restraint 
xis and CLHS is the Coefficient of 
ese same forces. In the printout 
e number of blades, which in this 
the solidity of the rotor system, 
blade stall indicator. It is 
ng the first azimuthal station 
80 degree and 360 degree position 
reater than the value of CQD at the 
muthal position. 
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8 . Harmon ic Anal ysis of 2 For ce 

The following is a harmonic analysis of the 
in the Z or shaft axis direction. Z force is posi 
the downward direction. This is PRINT Option 1000. 



f orces 
ive in 
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HARMONIC ANALYSIS OF DIMENSIONAL CZ 
A0+ A1 *COS (PSI)+B1*SIN(PSI)+A2*COS(2*PSI)+B2*SIN(2*PSI) 
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9 • Conv erged 31 a de Ana lysi s 

The next four pages contain the sample overall 
summary of events occuring on the blade at each azimuthal 
and radial position. This is PRINT Option 1. Variable 121, 
PPSI, controls the azimuthal intervals that are printed out. 



tput 


items 


are as follows. 


a . 


X - Center location of the blade segment 


b. 


ALPHA 


- Angle-of-Attack 


c. 


MACH 


- Local Mach number 


d. 


CL - 


Local Coefficient of Lift 


e. 


CD - 


Local Coefficient of Drag 


f. 


PHI - 


Local Inflow Angle 


g- 


L (LB/IN) - Lift produced per inch on segment 


h. 


Swee p 


- Sweep Angle of airflow 
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10. Harmonic of Air Loads 



The next printout contains the harmonic analysis of 
the lift generated per inch on each rotor blade segment. 
This is PRINT Option 10000. The airloads are computed in 
two harmonic forms, both containing terms out to and 
including the tenth harmonic. The forms are 

Lift per Inch = A0 - AlcosY - BlsinM? - A2cos2M> - 

B2sin2Y - A3cos3Y - B3sin3^ - AlOcoslOY - 

BIOsinlOY and 

Lift per Inch = A0 - C1sin(Y+4,) - C2sin(2Y+4 i ) 

-C3sin (3Y+<ty - ClOsin (10 y + ^) . 

The Ratio column is determined from the coefficients 
in Column C. From this ratio, one can immediately determine 
which airload harmonic is dominant and the relative 
relationship of this to the ether harmonic coefficients. 
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HARMONICS OF AIR LOADS 
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STATION 10 AO = 34.9860229 
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STATION 13 AO = 14.6387329 
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IV. GRP SAMPLE ANALYSIS 



The GRP program was executed using data representing a 
relatively new rotor blade. The results were compared with 
results predicted by the blade’s manufacturer. The rotor 
blade was an unsymmetrical blade. Ihe rotor radius was 
divided into three sections. Sections one and three were 
made of the same type airfoil. The blade included a sweep 
tip design. The blade and helicopter configuration analyzed 
are typical of a helicopter that could be used by the O.S. 
Navy in a LAMPS type mission. 

For the analysis, the blade was divided into 15 
segments. The program used the manufacturer's values for 
the First and Second Moment of Inertia about the Flapping 
Hinge, vice local blade mass densities. The GRP assumed 
uniform inflow for all flight velocities . It used a rigid 
blade analysis while the actual blade does have live twist. 
The program was run at five different flight weights, 
ranging from 16,359 to 20,829 pounds. A flat plate area of 
35.8 sguare feet was used at all speeds. The GRP was run 
for forward flight speeds of 40 to 160 knots at ten knot 
increments. The program was executed at sea level, tropical 
day condition. The manufacturer's predicted rotor 
horsepower was obtained from his Shaft Horsepower versus 
True Airspeed curves and was corrected to rotor horsepower 
by using the manufacturer's Mechanical Efficiency curves. 

The results of the analysis are shown in the next 
several tables. Table VI illustrates a comparison of the 
GRP required rotor horsepower divided by the manufacturer’s 
required rotor horsepower. The GRP agreed within an average 
of two percent on the entire range from 40 to 160 knots. 
The GRP agreed within an average of one percent for the 
cruise range between 70 and 140 knots. It can be seen that 
between 40 to 60 knots there is a much larger difference 
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between the two required horsepowers. It is felt that the 
inflow in this region is not uniform as assumed, but highly 
mixed and irregular. Also, the GRP results were less than 
the manufacturer's horsepower in the 150 to 160 knot range. 
This area represents the region of top speed for the 
helicopter, and much of the retreating blade is in the stall 
region. Also, it is expected that there is a change in 
fuselage attitude at this high speed, which would increase 
the flat plate area above what was used in the program. The 
GRP program's maximum endurance velocities agreed exactly 
with those predicted by the manufacturer. 





TABLE I 


WEIGHT = 16359 LBS 




VELOCITY 


GRP RHP 


MANUFACTURER'S RHP 


RATIO (GRP/MAN) 


40 


1105 


1152 


.96 


50 


1006 


1041 


.97 


60 


960 


973 


.99 


70 


955 


958 


1 . 00 


80 


984 


97 1 


1.01 


90 


1047 


1034 


1.01 


100 


1142 


1 109 


1 .03 


110 


1273 


1237 


1.03 


120 


1438 


1396 


1.03 


130 


1643 


1589 


1.03 


140 


1893 


1852 


1.02 


150 


2300 


2216 


1.04 


160 


2567 


2745 


. 94 





TABLE II 


WEIGHT = 17321 LBS 




VELOCITY 


GRP RHP 


MANUFACTURER* S RHP 


RATIO (GRP/MAN) 


40 


1185 


1247 


. 95 


50 


1075 


1122 


. 96 


60 


1018 


1038 


.98 


70 


1008 


1009 


1.00 


80 


1028 


1028 


1 . 00 


90 


1086 


1078 


1.01 


100 


1176 


1153 


1 .02 


110 


1303 


1264 


1 . 03 


120 


1465 


1445 


1.01 


130 


1669 


1633 


1.02 


140 


1921 


1902 


1 .01 


150 


2227 


2257 


. 99 


160 


2598 


279 1 


.93 




TABLE III 


WEIGHT = 19246 LBS 




VELOCITY 


GRP RHP 


MANUFACTURER'S RHP 


RATIO (GRP/MAN) 


40 


1366 


1449 


. 94 


50 


1229 


130 1 


. 94 


60 


1 1 48 


1186 


.97 


70 


1118 


1142 


. 98 


80 


1126 


1153 


. 98 


90 


1 172 


1199 


. 98 


100 


1255 


1273 


.99 


110 


1375 


1387 


. 99 


120 


1530 


1557 


. 98 


130 


1731 


1746 


.99 


140 


1983 


2038 


. 97 


150 


2290 


2440 


. 94 


160 


2665 


- 


- 
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TABLE 


IV 


WEIGHT = 19658 


i LBS 




iLOCITY 


GRP RHP 




MANUFACTURER' S 


RHP 


RATIO (GRP/MAN) 


40 


1407 




1496 




. 94 


50 


1264 




1340 




. 94 


60 


1178 




1231 




. 96 


70 


1142 




1173 




. 97 


80 


1149 




1 179 




. 98 


90 


1193 




1221 




. 98 


100 


1274 




1299 




.98 


110 


1392 




1414 




. 98 


120 


1546 




1566 




. 99 


130 


1747 




1773 




. 99 


140 


1999 




2048 




. 98 


150 


2305 




2477 




. 93 


160 


2682 




- 




- 




TABLE 


V 


WEIGHT = 20829 


LBS 




LOCITY 


GRP RHP 




MANUFACTURER' S 


RHP 


RATIO (GRP/MAN) 


40 


1532 




1627 




.94 


50 


1370 




1457 




. 94 


60 


1269 




1333 




. 95 


70 


1222 




1267 




. 96 


80 


1219 




1253 




. 97 


90 


1256 




1297 




. 97 


100 


1330 




1375 




. 97 


110 


1445 




1477 




. 98 


120 


1598 




1638 




. 98 


130 


1798 




1866 




.96 


140 


2045 




2157 




.95 


150 


2354 




2589 




.91 


160 


2747 




— 




- 
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TABLE VI 


RATIO COMPARISON 












WEIGHT 






RATIO 


VELOCITY 


16359 


17321 


19246 


19658 


20829 


AVERAGE 


40 


.96 


. 95 


.94 


.94 


. 94 


. 95 


50 


.97 


. 96 


.94 


.94 


.94 


.95 


.6 0 


.99 


.98 


.97 


.96 


. 95 


. 97 


70 


1.00 


1.00 


.98 


.97 


.96 


. 98 


80 


1.01 


1 . 00 


.98 


.97 


. 97 


. 99 


90 


1 .01 


1 .01 


.98 


.98 


.97 


.99 


100 


1 .03 


1.02 


.99 


.98 


. 97 


1 . 00 


110 


1 .03 


1.03 


.99 


.99 


.98 


1 . 00 


120 


1.03 


1 .01 


.98 


.99 


. 98 


1 . 00 


130 


1 .03 


1 . 02 


.99 


.99 


.96 


1 . 00 


140 


1 .02 


1 . 01 


.97 


.98 


.95 


. 99 


150 


1 . 04 


. 99 


.94 


.93 


.91 


. 96 


160 


.94 


. 93 


- 


- 


- 


.94 




AVERAGES FOR 


ENTIRE SPEED RANGE 






AVERAGE 


1 .00 


. 99 


.97 


.97 


. 96 


. 98 






AVERAGES 


FOR CRUISE 


RANGE 










70 


- 140 KNOTS 








AVERAGE 


1 .02 


1 . 01 


.98 


.98 


. 97 


. 99 
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V. CONCLUSIONS 



The logic and theory used in the GRP program was 
investigated and found to be sound. However, there were three 
discrepancies in the Navy's version of the program that did 
require attention. The calculations in the reverse flow section 
on the rotor were incorrect, the calculation of the chord at the 
75 percent radius position was incorrect and the original Trim 
Option for reducing moments would not work. All of the above 
discrepancies were corrected. 

Three desirable features were added to the GRP program. 

First, the ability to analyze a rotor blade composed of more than 
one airfoil type was added. The program will now accept up to 
five different airfoil data input decks for use in analyzing a 
rotor system. Secondly, the program would only calculate performance 
in level flight. The ability to calculate performance in climbs 
and descents has been added. Lastly, the program will now calculate 
the aerodynamics for a swept tip rotor blade design. 

The results of the sample analysis described in Section IV 
indicates that the program does produce highly accurate performance 
predictions. The averaged GRP rotor horsepower was within two 
percent of the manufacturer's data. The results were within one 
percent when compared in the area of normal cruise flight. The 
GRP, in this analysis, assumed uniform inflow, constant flat plate 
area and a rigid rotor blade. It is felt that while complicated, 
computer- time-comsuming procedures can be taken to reduce these 
assumptions, they are not warranted if the GRP is to be used 
strictly as a helicopter performance prediction program. 
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